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ABSTRACT 

An  empirical  study  of  the  time  and  space  variability  of  the  wind  as  measured  by  a  tower-mounted  wind 
instrument  and  triple-theodolite  pilot  balloon  observations  at  Green  River,  Utah,  during  the  months  of 
June  and  July  1964,  is  presented. 

Simultaneous  wind  measurements  at  500  ft  above  the  ground  are  rNmpared.  Mean  differences  in  wind 
direction  was  15.0  deg  with  a  sample  standard  deviation  of  12.2  deg.  The  mean  speed  difference  was  6.2  ft 
sec-1  with  a  sample  standard  deviation  of  5.5  ft  sec-1. 


1.  Introduction 

Wind  measuring  instruments  have  been  mounted  on 
towers  to  vary  ing  altitudes  up  to  approximately  1500  ft 
in  order  to  obtain  a  wind  profile.  Another  system  usetl  to 
obtain  winds  in  this  layer  is  the  pilot  balloon  observa¬ 
tion  technique  where  one  or  more  theodolites  track  a 
balloon.  Since  the  first  sy  stem  is  in  essence  an  Eulerian 
type  measuring  system  and  the  balloon  system  is  a 
semi-Lagrangian  measurement,  the  question  arises  as 
to  the  compatibility  of  the  wind  profiles  generated  by 
these  systems.  That  is,  if  both  systems  are  measuring 
the  wind  perfectly,  then  the  wind  should  be  exactly 
the  same  at  any  given  point  along  the  profile. 

Experimentally  it  has  been  found  that  large  differ¬ 
ences  are  present  in  the  profiles  generated  bv  the  two 
systems  involved.  It  is  realized  that  many  factors  can 
influence  the  resultant  wind,  among  which  are  space  and 
time  variability,  terrain  effects,  systems  error,  balloon 
aerodynamics,  tower  shadow  effect,  and  reduction 
processes.  It  is  beyond  the  scope  of  this  paper  to  at¬ 
tempt  to  determine  the  extent  of  influence  of  any  one  of 
these  factors  on  the  wind  profile.  The  purpose  of  this 
study  is  to  show  the  “apparent”  wind  differences  at  a 
given  level  of  the  profile  without  attempting  to  deter¬ 
mine  the  cause.  The  level  selected  for  this  study  was  the 
5<J0-ft  level  utilizing  a  three-cup  anemometer  and  a 
three-theodolite  wind  measuring  system. 

2.  Data  collection  and  reduction 

The  500-ft  tower  is  located  approximately  400  ft 
from  a  missile  launch  pad  at  Green  River,  Utah.  Only 
the  data  from  the  anemometer  at  the  500-ft  level  were 
utilized  in  this  study  since  the  problem  concerns  the 
differences  in  the  measured  wind  at  the  top  of  the  tower 
(an  Eulerian  system)  and  the  pibal  measured  wind 
(approximate  Lagrangian  system)  at  the  same  level. 
The  data  from  the  anemometer  were  transmitted 


electronically  to  a  digital  computer  at  a  sampling  rate 
of  one  point  per  second  and  averaged  over  a  10-sec 
period  centered  around  the  time  at  which  the  pilot 
balloon  reached  a  height  of  500  ft  to  provide  a  mien 
wind  for  comparison.  The  thickness  of  the  layer  trav¬ 
ersed  in  a  10-sec  period  varies  with  ascent  rate  but  is 
generally  140  to  180  ft. 

The  pibal  system  consisted  of  three  semiautomatic 
theodolites  placed  as  an  equilateral  triangle  with  sides 
of  approximately  5345  ft.  In  this  system  the  balloon 
was  visually  tracked  by  an  observer,  and  the  azimuth 
and  elevation  angles  were  electronically  transmitted 
on  a  real  time  basis  (Duncan,  1963)  to  a  digital  computer 
at  White  Sands  Missile  Range,  providing  a  wind  value 
each  second.  These  w  ind  data  were  also  averaged  over 
a  10-sec  period  centered  around  the  time  the  balloon 
passed  through  the  500-ft  level.  The  pibal  winds  were 
calculated  utilizing  Duncan’s  mathematical  derivation 
which  is  presented  in  the  Appendix. 

3.  Discussion 

The  tower  wind  is  a  measurement  at  a  fixed  point. 
This  wind  is  averaged  over  a  10-sec  interval.  The  pibal 
wind  is  a  mean  wind  also  over  a  10-sec  interval  (a  layer 
approximately  170  ft  thick)  centered  at  the  500-ft 
height  which  involves  horizontal  as  well  as  vertical 
space.  It  has  been  stated  (Panofsky  and  Lumley,  1964) 
that  Lagrangian  and  Eulerian  wind  variances  are  theo¬ 
retically  equal  in  stationary,  homogeneous,  incompres¬ 
sible  turbulence.  A  comparison  between  tower  and  air¬ 
plane  wind  measurements  at  heights  of  295  ft  and  394  ft 
by  Lappe  el  al.  (1959)  showed  that  the  theory  appears 
to  be  valid  for  horizontal  wind  fluctuations. 

The  surface  roughness  at  Green  River  (Armendariz, 
1965)  precludes  homogeneity  of  the  atmosphere  in  the 
lower  layer.  It  is  felt  that  considerable  mechanical 
turbulence  exists  with  a  dominance  of  eddy  sizes  of  a 
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Fjc.  1.  A  typical  low  level  wind  profile  comparison  lietween  pibal 
and  tower  wind  measurements  at  Green  River,  Utah. 


scale  less  than  the  distance  between  the  tower  and  the 
balloon. 

In  a  study  by  Moses  and  Daubek  (1961)  it  was  con¬ 
cluded  that  towers  distorted  the  wind  flow  thereby 
producing  errors  in  wind  velocity  measurements  ob¬ 
tained  by  anemometers  mounted  on  them.  In  some  cases 
they  found  that  the  tower-mounted  anemometer  gave 
wind  speed  readings  appreciably  higher  than  a  reference 
anemometer,  and  in  other  cases  for  given  wind  direc¬ 
tions,  there  was  a  substantial  reduction  in  the  wind 
speed  recorded  from  the  tower  as  compared  to  the 
reference  anemometer.  This  tower  wind-shadow  effect 
reduced  the  wind  measurements  nearly  one  half  in 
light  winds  and  nearly  25  per  cent  for  speeds  between 
10—15  mph.  Preliminary  results  in  the  wind  tunnel  at 
Colorado  State  University  (Cermak,  1964)  showed  a 
reduction  of  10-15  per  cent  for  the  same  speed  range. 
There  is  also  some  evidence  that  the  ordinary  pilot 
balloon  does  follow  the  actual  wind  flow  (MacCready 
cl  at.,  1964;  Scoggins,  1965;  Killen,  1960)  and  therefore 
is  not  a  true  Lagrangian  system. 

4.  Results 

The  mean  direction  difference  (Table  1)  for  the  71 
cases  was  15.0  deg,  with  a  sample  standard  deviation 
of  12.1  deg.  The  greatest  mean  difference  recorded  was 
on  12  June  when  the  mean  difference  was  25.5  deg 
and  the  least  was  on  14  June,  a  mean  of  7.1  deg. 
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Fkj.  2  Correlation  of  tower  wind  direction  at  500  ft  and  pibal  wind  direction  at  the  same  level  at  Green  River,  Utah, 

during  June  and  July  1964, 
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Fio.  5.  Relation  between  the  difference  in  wind  speed  of  tower 
and  pibal  at  500  ft  to  distance  from  balloon  to  500-ft  anemometer. 
Green  River,  Utah  71  cases. 


The  mean  wind  speed  difference  was  6.2  ft  sec-1,  with 
a  sample  standard  deviation  of  5.5  ft  sec-1.  The  greatest 
mean  speed  difference  occurred  on  17  June  with  a  mean 
speed  difference  of  14.1,  whereas  on  12  June  the  mean 
was  only  2.1  ft  sec-1. 

The  linear  correlation  of  the  speed  and  direction 
between  the  tower  and  the  pibal  wind  was  r„  =  0.86  and 
rp=0.95,  respectively. 

It  is  interesting  to  note  for  the  case  illustrated  in 
Fig.  1  that  the  direction  of  the  east-west  component 
for  the  pibal  is  generally  opposite  to  that  of  the  tower, 


Table  1.  Pibal  and  tower  wind  comparison  at  500  ft.  Green 
River,  Utah.  Direction  difference  is  in  deg;  speed  difference  in 
ft  see-1. 


Dates 

Observa¬ 

tions 

Wind  direction 
Dr-Dp  SiOf-Df) 

Wind  speed 

Ut~  Up  S(Vt-v 

10  June  1964 

10 

9.2 

5.7 

ii.i 

5.0 

12  June  1964 

12 

25.3 

16.4 

2.1 

1.9 

14  June  1964 

12 

7.1 

6.9 

3.0 

2.6 

15  June  1964 

11 

17  9 

12.4 

6.7 

3.0 

16  June  196# 

9 

19.6 

11.6 

24 

2.7 

17  June  1965 

12 

11.0 

5.3 

14.1 

4.9 

8  July  1964 

5 

18.4 

4.6 

2.4 

2.2 

Total 

71 

{Dr ‘-Dp)  ti  —  15.0 

5er”6.1  ft  sec-1 

5crj>“10-4  ft  »*c_' 

(VT-V,)n-  62 

Soy  — 40.4° 

5.5 

So, -48.5° 

0r-14.6ft.ee-1 

r,—  +0.86 

f?r- 20.2  ft  tec-1 

fp"  +0.95 

0r-2O5° 

S» sample  standard  deviation 

/>r-213* 

Dr-wind  direction  of  tower 
Dr— wind  direction  of  pibal 

U- wind  speed 
r- correlation  function 

end  the  speed  of  the  north-south  component  is  less  for 
the  tower  than  for  the  pibal.  These  differences  are  indi¬ 
cated  in  Figs.  2  and  3. 

Fig.  2  shows  the  majority  of  the  wind  directions  plot 
to  the  right  of  the  1:1  correspondence  line.  In  Fig.  3 
most  of  the  wind  speeds  are  also  to  the  right  of  the  1:1 
correspondence  line,  indicating  that  the  tower  is 
generally  reporting  lower  wind  speeds. 

Fig  4  is  a  plot  of  the  tower  wind  speed  as  abscissa 
and  the  ratio  of  the  tower  wind  to  pibal  wind  as  ordinate. 
In  59  of  the  71  cases  the  ratio  was  less  than  one,  indi¬ 
cating  a  lower  value  for  the  reported  tower  wind  than 
for  the  pibal.  The  ratio  is  significantly  less  than  unity 
for  southeasterly  wind  speeds  of  less  than  15  ft  sec'1, 
and  for  southwesterly  winds  with  speeds  greater  than 
15  ft  sec-1. 

Fig.  5  is  a  scatter  diagram  of  the  differences  between 
the  tower  and  the  pibal  wind  speeds  as  a  function  of  the 
balloon  distance  relative  to  the  tower  at  the  time  the 
balloon  reached  the  500-ft  level.  As  would  be  expected, 
an  increased  wind  speed  difference  with  distance  is 
indicated. 

5.  Conclusion 

Differences  exist  in  wind  observations  obtained  with 
tower-mounted  anemometers  and  with  pilot  balloons 
at  the  500-ft  level.  Mean  direction  difference  was  15.0 
deg  and  the  mean  speed  difference  was  6.2  ft  see1. 
Good  agreement  between  the  wind  directions  from  the 
tower  and  pibal  measurements  is  evident  from  a  linear 
correlation  of  0.95.  In  59  of  the  71  cases  investigated  at 
Green  River,  Utah,  the  tower  recorded  a  wind  speed 
lower  than  that  obtained  from  the  pibal  observation. 
It  is  believed  that  these  differences  can  be  attributed 
to  some  combination  of  the  following  factors:  space 
and  time  variability;  tower  shadow  effect;  balloon 
aerodynamics;  terrain  effects;  systems  error;  and 
reduction  process. 


APPENDIX 

A.  Determination  of  Balloon  Position 

The  position  of  a  point  in  space  ran  be  determined 
by  an  azimuth  and  at,  elevation  angle  from  one  known 
location  and  an  azimuth  or  elevation  angle  from  another 
known  location.  However,  since  no  tracking  system  is 
error-free,  a  more  nearly  correct  solution  can  be  ob¬ 
tained  by  increasing  the  amount  of  angular  data. 
Statistical  procedures  are  then  applied  to  these  data 
to  obtain  the  most  probable  position. 

The  techniques  used  in  this  system  determine  the 
point,  assumed  to  be  the  balloon  position,  such  that 
the  sum  of  the  squares  of  distances  from  this  point  to 
each  line  of  sight  (as  determined  by  the  azimuth  and 
elevation  angle  from  the  theodolite)  is  a  minimum. 

The  position  is  calculated  with  respect  to  an  earth- 
fixed  Cartesian  coordinate  system,  OXYZ,  defined  as 
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The  function  to  be  minimized  is 

F~Zh,*.  (4) 

<-i 

By  using  the  formula  for  the  distance  between  two 
points  and  the  properties  of  direction  cosines,  we  have, 
after  simplification 

f-z  Kx-x.Y+(y-riY+b-z.y 

-[«i(r-.Y,H-61fy-ri)+c<(*-/.-)J}.  (5) 
A  necessary  condition  for  F  to  be  a  minimum  is 
d  F  BF  BF 

- :  —  =  —  =  0.  (6) 

Bx  By  Sz 

Performing  the  operations  in  Eq.  6,  we  have,  after 
considerable  simplification 

A  ii*+.l  uy~i~  A  iiz= Di, 

Anx+Any+A  (7) 

A  n.r + A  ny+-l  i&  =  Di, 

where 

.1„  =  £  (1  — a,1) 


follows :  Z  lies  along  the  zenith,  X  and  Y  are  oriented 
such  that  OXYZ  is  a  tight-hand  system  (see  Fig.  6). 
In  practice  either  X  or  Y  is  usually  oriented  positive 
northward ;  however,  such  a  restriction  is  irrelevant  to 
the  solution. 

The  following  notations  are  defined : 

P,=  (X is  the  position  of  the  >th  theodolite. 
(.r,y,z)  is  the  point  to  be  determined, 
o,  is  the  azimuth  angle  measured  positive  clockwise 
from  a  line  through  Pi  parallel  to  the  A'  axis  to  the 
projection  on  the  A--!'  plane  of  the  line  of  sight  from 
the  ith  theodolite. 

6,  is  the  elevation  angle  from  the  ith  theodolite 
measured  positive  upward. 
d,  is  the  distance  from  (.r,y,z,l  to  (XitYi,Zi). 
h,  is  the  distance  from  (x,y,z,)  to  the  line  of  sight 
from  the  ith  theodolite. 

t,  is  the  angle  between  d,  and  the  ith  line  of  sight, 
f-.li.f.  are  the  direction  cosines  of  d,. 
a„b„c,  are  the  direction  cosines  of  the  ith  line  of  sight. 
From  Fig.  6  it  is  clear  that 

x—  X,  y—Y,  s-Zi 

- ,  e.— f.-~— ;  (1) 

di  d%  d% 

a4»cosfl(Cosa„  fc.^cosfl,  sina,,  r(  =  sin04,  (2) 

and 

A,=(/iSin«(.  (3) 


iu=E  (i-*.1) 

.-I 

A»=Z  (1-c*) 

i-l 

n 

*1 1*“  A  t\  —  —  £  atbi 
•— » 

N 

.1  1.1=  .1  31  =  —  £  0* f * 
i-l 

*  1 23  =  A  as  =  —  £  , 

»— i 

Di  =  2[m(l  —  a,’)Xi  —  a,b,)',—  a,cyZ,2 

nt=Z£(-a,b,X,+  (l-b,')Yi-b.<-J,] 

D,=Sl-airiX,~biCiYi+  (1-c  *)/,]. 

The  solution  of  the  system  of  Eqs.  (7)  is  the  desired 
point  (r,y,z,V 

B.  Calculation  of  the  Wind  Components 

The  wind  profile  is  determined  by  numerical  differen¬ 
tiation  of  the  position  profile;  the  differentiation  is 
performed  componentwise.  In  general,  the  particular 
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numerical  differentiation  technique  one  uses  depends 
on  several  fac  tors  which  are  peculiar  to  his  particular 
problem.  After  consideration  of  data  available,  the 
accuracy  obtainable,  and  the  amount  of  computation 
required,  we  have  chosen  the  following  procedure. 

Let  A  be  an  altitude  at  which  the  wind  calculation  is 
desired.  Let  (Ao,yo,2o)  be  the  point  on  the  position  pro¬ 
file  such  that  is  nearest  to  h.  Eleven  successive  (time- 
wise)  position  points  are  chosen  with  (jro,yo,2o)  as  their 
temporal  midpoint.  For  notation  purposes  these  points 
are  given  subscripts  (—5,  —4,  .  .  .0,  .  .  5)  where 

increasing  subscript  indicates  increasing  time.  A  parab¬ 
ola  is  fitted  to  these  points,  componentwise,  and  is 
differentiated  to  obtain  the  balloon  speed,  the  negative 
of  which  is  assumed  to  be  the  wind  speed  (the  direction 
from  which  wind  is  blowing). 

The  mathematics  involved  is  a  straightforward  appli¬ 
cation  of  least-squares  techniques.  Let  To  be  the  time 
associated  with  the  midpoint,  and 

t-t* 

t  =  — - — . 

A/ 

It  suffices  to  find  the  coefficients  of 


x=a<,+  ait+aiP. 


(8) 


This  is  considerably  simplified  since  it  is  required  to 
evaluate  only  the  derivative  at  T-To,  i.e., 


d  a  i  a, 

x  =  —  (a0+<i  d +ajff)  =  — F  2 — /, 
dT  At  At 


and 


x|To=- 


ai 

M 


(9) 

(10) 


One  of  the  equations,  readily  obtainable  by  classical 
least-  squares  techniques,  is 


(ID 


We  observe  that  Zt  =  ZP  =  0.  Hence, 


Xtx 
ai  =  — . 

ZP 

Thus 

a  i  Zlx  1  s 

r  =  —  = - = -  51  »■*'. 

At  A tZP  110A/--S 


(12) 

(13) 


is  the  .r  component  of  the  balloon  speed.  Thus, 


and 


1  5 

H  ,  = -  E  ix„ 

llOAt  i— ■> 

1  » 

trv= -  E  ty,. 

llOAt  6 


(14) 
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